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Abstract — Decode-and-forward (DF) cooperative communica- 
tion based on free space optical (FSO) links is studied in this 
letter. We analyze performance of the DF protocol in the FSO 
links following the Gamma-Gamma distribution. The cumulative 
distribution function (CDF) and probability density function 
(PDF) of a random variable containing mixture of the Gamma- 
Gamma and Gaussian random variables is derived. By using the 
derived CDF and PDF, average bit error rate of the DF relaying 
is obtained. 

Index Terms: Bit error rate, decode-and-forward relaying, 
free space optical links, Gamma-Gamma fading. 

I. Introduction 

Free space optical communication has recently drawn sig- 
nificant attention because of very large bandwidth, low imple- 
mentation cost, and excellent security 0, 0. The intensity 
modulation and direct detection (IM/DD) with on/off keying 
(OOK) are used for simplicity in the FSO communication. 
However, in the OOK system it is difficult to choose an 
optimal detection threshold in presence of the atmospheric 
turbulence . In order to enhance the performance of the FS O 
system, the subcarrier intensity modulation (SIM) is proposed 
in 0, 0. The SIM leverages on advances made in signal 
processing as well as the maturity of radio frequency (RF) 
devices such as highly selective filters and stable oscillators, 
which then permits the use of modulation techniques such as 
phase-shift keying (PSK) and quadrature amplitude modula- 
tion (QAM) 0. 

Cooperative communication is a technique to improve the 
diversity and coverage area in a distributed antenna based 
RF communication system For FSO systems, coop- 

erative communication has been studied in ||8l- lfT0l , where 
the turbulence-induced fading in the FSO links is modeled 
by the log-normal distribution which is commonly used to 
model weak turbulence conditions. However, the Gamma- 
Gamma distribution has been widely adopted for study of 
behavior of the FSO links under wide range of atmospheric 
turbulence conditions (weak to strong) because it fits well to 
the experimental results 0-0, 1(171. ifTII. 

In this letter, we analyze the average bit error rate (BER) 
performance of the DF based cooperative FSO communication 
system under the Gamma-Gamma fading channels. We adopt 
the SIM scheme with the binary phase shift keying (BPSK) 
modulation. We first derive the CDF and PDF of a random 
variable (RV) containing mixture of the Gamma-Gamma and 
Gaussian RVs. Then, with help of the CDF and PDF expres- 
sions, we derive a power series expression of the average BER 
of the DF scheme. 
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II. System Model 

Let us consider a subcarrier intensity-modulated FSO co- 
operative system containing a source terminal (S) with two 
transmit apertures (for transmission of the data to the relay 
and the destination), one destination node (D) with two receive 
apertures (for receiving the signals transmitted by the relay 
and the source), and a single relay (R) with one transmit 
aperture (for transmitting signals to the destination) and one 
receive aperture (for receiving signals from the source). The 
transmission from S to D is performed in two phases. In first 
phase, S transmits BPSK signal x € {1, —1} to R and D by 
using different transmit apertures. The relay demodulates the 
data of S in symbol-wise manner. It is assumed that S uses 
cyclic redundancy check (CRC) code such that R can know 
whether it has demodulated the data correctly or not. The relay 
transmits the demodulated symbols to D in next phase only if 
it has decoded the data of S correctly. The data received in R 
and D during the first phase will be 

y s ,d = Vs,dI s ,dX + e S} d, 

where I s r with E {I SJ -} = 1 and I s _d with E {I s d} = 1 
are the real-valued irradiance of the S-R and S-D links, 
respectively, following the Gamma-Gamma distribution, e S;r 
and e s ,d are signal-independent complex-valued additive white 
Gaussian noise (AWGN) with zero mean and variances N s , r 
and N Si d, respectively, and r] s d and r/ s . r are the optical-to- 
electrical conversion coefficients. If R decodes the data of the 
source correctly, then data received in D in the next phase 
from R will be , 

Ur,d = Vr,d-lr,dX + e r ,d, (2) 

where I r ,d with E {I r ,d} = 1 is the real-valued irradiance from 
R to D, following the Gamma-Gamma distribution, e r ^ is the 
signal-independent complex-valued AWGN noise with zero 
mean and N r> d variance, and r] r d is the optical-to-electrical 
conversion coefficient. 

The destination uses the following decoder: 

U + vV ^ 0, (3) 
-l 

where U 4Re jy* dVs^d/^d}' V 

4Rs ^y* d rj r ^dI r ,d/N ri d^, v = 1 if R transmits, and v = if 
R remains silent. 

III. Characterization of A Random Variable 
containing mixture of gamma-gamma and 
Gaussian Random Variables 

Let us define a random variable Y = aZ 2 + bZE, where a 
and b > are constants, and Z is Gamma-Gamma distributed 

as 
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where K c (-) is the modified Bessel function of the second 
kind of order c (131, the parameters a and j3 are related to the 
atmospheric turbulence conditions ||3|, |[T2"1 . and E ~ Af (0, 1) 
is the real-valued Gaussian distributed RV. 

Let W = Z 2 , then the PDF of W can be easily obtained 
from as 
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By using the series expansion of K a _p (2 v / a/3w 1 / 4 ) Ifl3l 
Eq. (9.6.2) and (9.6.10)] in © we have 
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Theorem 1: The CDF ofY is given as 
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anfif B (■,■) is the Beta function HT4\ Eq. (8.380.1)]. 
Theorem 2: The PDF of Y is given as 
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Refer Appendix J] for proof of Theorem 1 and 2. 



IV. Bit Error Rate of DF Relaying 

Let x = 1 be the transmitted bit by S but R decides that 
— 1 was transmitted, then the probability of error of the S-R 
link can be written as 

P e S ' R = Pr (X < 0) , (18) 
where X — 4Re {y* r r) s ,rIs,r/Ns, r } has the same distribution 
as that of Y with a = a s>r = Ar] 2 r /N s _ r and b = b s r = 
2\/2r) sr / y/N s-r . Therefore, from (|8) and (TT8l , we have 

P!' R = I (a 3tr , /3 s , r , a 8 , r , b„, r , 0) , (19) 
where {a Sjr ,/3 s r } are the turbulence parameters of the S-R 
link. The probability of error in wrongly decoding 1 as -1 in 
the destination receiver can be written by using ((3} as 

P e D = P e SR Pr (U < 0) + (1 - P e SR )Pr (U + V<0), (20) 
where the RVs U and V have the same distribution as that of Y 
with I a = a M = 4rj 2 jN Stdl b = b %d = 2y/2r] s4 / y/T^j 

and ja = a r4 = 4r]^ d /N r , d , b = b r4 = 2\/2<7r,d/ y/N rjd j, 
respectively. We can rewrite (l20t as 

pD = pS.RpS.D + £ _ pS,R)p r (f/ + < 0) , (21) 
where P| ' D can be obtained similar to ( fT9l and it is shown in 
Appendix HJ that for a s ,db 2 d = a r ,db 2 d , we have 
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Pi(U + V < 0) = ^ (Pr,d,a r ,d)d k2 V3 a , d ,a a , d ) 
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xd k2 {a s ,d,P s ,d)M(g,Pr,d/2 + k 1 /2,a s . d /2 + k 2 /2)}, (22) 
where g = [a r , d , b r , d , a s>d , b Sid ], and {a r , rf , /3 rjrf } and 



1} are the turbulence parameters of the R-D and S-D 



links, respectively; 
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where g = [ai, 61, a 2 , 6 2 ], 2 Pi (•,-;•;•) is the Hypergeometric 
function JH), and d kl (a r ,d, Pr,d) and dk 2 {a s ,d, Ps,d) are 
given in (0. 

V. Numerical Results 

The analysis and simulation is done under strong (a = 4.2, 
P = 1.4) and moderate (a = 4.0, P = 1.9) atmospheric 
turbulence conditions |@J. It is assumed that S-R, S-D, and R- 
D link have the same average signal-to-noise ratio (SNR) and 
atmospheric turbulence conditions and r) s , d = i] s .r = Vr,d = L 

In Fig. [U we have plotted BER versus SNR performance 
of the three node cooperative FSO system employing SIM 
with BPSK and using DF protocol with selective relay which 
transmits only if it decodes the data of S correctly and 
perfect relay which commits no error in detection, in strong 
and moderate Gamma-Gamma fading channels which vary 
independently in every time-interval. The SNR on y-axis of 
Fig.[T]indicates the average SNR of the S-D link. We have also 
plotted BER versus SNR performance of a non-cooperative 
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Fig. 1 

Average BER of the DF relaying with SIM BPSK in 
Gamma-Gamma fading channels. 

direct transmission based system, where S transmits at double 
power than the cooperative FSO system. It can be seen 
from Fig. Q] that the cooperative FSO system significantly 
outperforms the direct transmission based system under strong 
and moderate turbulence conditions. For example, the selective 
relay based DF FSO cooperative system can achieve an 
SNR gain of 14 dB over the non-cooperative scheme for 
a moderate turbulence regime at BER=10~ 4 . Moreover, the 
selective transmission based DF scheme have much better 
diversity than the direct transmission scheme and works very 
close to the DF scheme with perfect relay under strong and 
moderate turbulence conditions. 

VI. Conclusions 

We have analyzed the DF based FSO cooperative system 
over the Gamma-Gamma fading channels. The advantages of 
using cooperative communication over a direct transmission 
based non-cooperative system are established by using analysis 
and simulation. It is shown that the BER of the DF FSO system 
significantly reduces at high SNR as compared to the non- 
cooperative system. 

Appendix I 
Derivation of CDF and PDF of Y 

Since f Y \Zab{y) ~ M (aW,b 2 W), the CDF of Y, i.e., 
F Y (y) = Pr(Y<y) = E w {1 - Pr (Y > y\W, a, &)}, can 
be written by using 03] Eqs. (2.3.12) and (2.3.14)] as 

where Q (•) is the q-function 03] Eq. (2.3.10)]. From ©, 0, 
and d24"i i. we get ^ in terms of / (a, (3, a, b, y) defined in (O, 
where 
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By using the following relation: Q (x) 



dw, to > 0. (25) 



where erfc(-) is the complimentary error function, and JT6 
Eq. (2.8.4.3)] in (J25J, we get dTUJl and O- By subtitut- 
ing y = and using the following relation: Q (x) = 



1 rir/2 ,„ 
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Eq. (3.191.3)] in (J25J, we get {□}. 



> 0, 01 Eq. (3.381.4)], and 01 



By taking partial derivative of (|8) with respect to y, we 
get (13[ in terms of B (a, (3, a, b.y) defined in ( [Ti l, where 



D (a, b, to, y) = — J (a, b, m, y) 
by 



(26) 



The right hand side of d26l l can be simplified by us- 
ing 423, 01 Eq. (0.410)], and the relation ^Q(x) = 
— -JL= e - x ~/ 2 . Then we can obtain ( [131 ) and ( fTTI i by using 01 
Eq. (3.471.9)], and dT6b can be obtained by substituting y = 

and using 01 Eq. (3.381.4)]. 

Appendix II 

Derivation of Pr (U + V < 0) 
After some algebra, we can write 

° /-oo poo 

Pr(U + V < Q)= Fu(v)f v {-v)dv+ Fu(-v)f v (v)dv. (27) 
Jo Jo 

Let us define the following function: 
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K mi -i/ 2 (aiy/bf) K m2+1/2 (a 2 y/bl) dy. (28) 



Equation (l22l i can be obtained by using (0, ( fT3l . (l27T i. ( l28l . 
and QU Eq. (2.16.33.1)]. 
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